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Abstract 
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Tlj" ' The differential equation governing distance-redshift for partially filled-beam optics in pressure- 

(N : 

free Friedmann-Lemaitre-Robertson- Walker Cosmology (FLRW) is shown to be the Lame' equa- 
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tion. The distance-redshift, D{z), discussed is appropriate for observations in inhomogeneous 
_ cosmologies for which lensing by masses external to the observing beam are negligible and for 

00 

f-*) ■ which lensing by transparent matter within the beam can be approximated by a homogeneous 

m '. 

■ mass density expanding with the FLRW background. Some solutions of the derived Lame' equa- 



tion are given in terms of Weierstrass elliptic integrals. A new simplified and useful expression for 



O ' filled-beam D{z) in standard flat FLRW is also given. 
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I. INTRODUCTION TO GRAVITATIONAL OPTICS 



One constant goal of cosmology has been to determine the large scale geometrical struc- 
ture of the Universe, i.e., the Hubble parameter Hq and the mass density parameter Q^. 
More recently an additional dynamical parameter Q\ (a cosmological constant or vacuum en- 
ergy term) has been added to the "must know" list and even more exotic forms of matter and 

n 

associated parameters are being speculated about, e.g., quintessence, see j30||. One statisti- 
cally significant constraint on these parameters comes from fitting the observed magnitude- 
redshift of type la supernovae to Hubble curves predicted by the FLRW models. 
Magnitude- redshift m{z), or equivalent distance- redshift D{z), for FRW (same as FLRW 
but without A) was first given by Mattig ^ . In contrast to the real Universe, FLRW mod- 
els are homogeneous and isotropic on all scales. The first modification to the theoretical 



50|. This 



Hubble curves to account for the lack of homogeneity dates back to Zel'dovich 
modification amounted to correcting for the decreased converging effect of the decreased 
mass density in the line of sight between the observer and the distant object. The idea for 
this modification is that: (i) inhomogeneous matter exists in the universe at the expense 
of homogeneous matter, pj = po — pn, (h) the effects of inhomogeneous mass should be 
accounted for by gravitational lensing. It is this same correction to the Hubble curve being 
discussed here. 

The original geometrical construction by Zel'dovich was for a flat A = universe in which 
all matter was absent from the observing beam, pn = 0. He was the first to give a differential 
equation from which a corrected distance-redshift was derived. The cosmological constant. A, 
was noticeably absent in most calculations done in those days; it was considered to be nothing 
more than a mistake made by Einstein. The first derivation of the differential equation for 
D{z) with 7^ 1 was by Dashevskii and Zel'dovich [8J but again only for the empty beam 
case. It was Dashevskii and Slysh [3| that extended the equation to partially filled observing 
beams, p/f 7^ 0. The differential equation evolved as the number of parameters increased 
and the types of cosmologies changed, see Kantowski and for some history of solutions 
see I23 I. 

n 

The first attempt to make use of the corrected magnitude-redshift by Bertotti [3J was 
met with widespread skepticism. The primary objection (which is valid for weak lensing 
but not for strong) was that the mean magnitude of a distant object didn't agree with the 
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standard homogeneous results of Mattig [31[ . The primary source of the 'error' was thought 
to be the perturbative nature of the cosmological model used. In Betrotti [3| part of the 
homogeneous FRW mass density was removed and the effect of this missing mass density 
on optical observations was accounted for in accordance with the earlier work of Zel'dovich 
[soj l . Additional point masses were superimposed to account for lensing by nearest neighbor 
galaxy impacts. Another more comprehensive numerical calculation by Gunn did not 
contradict this previous work. 

To overcome the suspected approximation error, the author was encouraged to develop 
geometrical optics for a type of universe called Swiss cheese. This universe is homogeneous 
FLRW with numerous spherical holes removed and replaced by point masses Schwarzschild 
solutions or other spherically symmetric matter |4], (see [28] for a history of these 
solutions) appropriately constrained to satisfy Einstein's gravity theory. Because this model 
was an exact solution to Einstein's GR, questions related to perturbative errors disappeared. 
With the introduction of the complex Swiss cheese model came the necessity of following 
waves as they pass from homogeneous FLRW into and through inhomogeneous Schwarzschild 
or Bondi and out again and so on. With Swiss cheese it could be demonstrated that a 
wave going through an inhomogeneity re-emerged with an insignificantly modified redshift- 
radius relation {1 + z = Rq/R) and that the affine parameter value was likewise unaltered 
from FRW. These results were necessary, but were suspect assumptions made when using 
superposition models. Geometrical optics in Swiss cheese could be done only because of the 
earher work of Sachs 0] and Jordan, Ehlers, and Sachs [2^ that developed the optical-scalar 
equations for light propagation in any gravity field. 

The geometrical optics approximation for light waves is equivalent to a bundle of non- 
rotating rays traveling from a source in the distant past to an observer here and now. It is 
the fractional increase in the area A of this beam that determines the fractional decrease in 
the luminosity of the object as seen by us, i oc 1/A. What complicates evaluation of the area 
A of an observing beam is the existence of two accompanying properties, the distortion of 
the beam (measured by its shear a or equivalently = cr A used below) and the orientation 
of the distortion (measured by an angle (p). The dynamics of these three parameters are 
intricately coupled as can be seen in the following form Kantowski |2l| of the optical scalar 
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equations of 



39| and 




^' = AC cos{p - (p) 



^'^ = AC sm{p -(f)). 



(1) 



These equations are complicated in their details but are easy to understand in principle. 
They contain three independent variables A, ^, and (p that represent properties of a particular 
geometrical optics wave. The other three terms TZ, C, and f3 are driving terms which are 
properties of the gravitational field in which the wave propagates. To use these equations 
first connect the source and observer with a single central light ray, a null geodesic of the 
gravity field, parameterized by an affine parameter, i.e., the prime derivatives in (PJ. You 
then start the three independent wave variables at the source with initial data appropriate 
for a point source in Euclidian geometry, A oc r^, ^ = 0, and = arbitrary. As the bundle 
of rays propagates through the gravity field they are converged {A is decreased) by the Ricci 
tensor (the negative TZ term) and are distorted or sheared (i.e., ^ = a A becomes non-zero) by 
the Weyl tensor C (the conformal curvature part of the Riemann tensor). The Ricci tensor is 
determined by the local density of matter within the bundle of rays whereas the Weyl tensor 
is determined by inhomogeneous matter near but external to the rays. The orientation of 
the distortion of the rays (p is driven by the orientation j3 of the inhomogemeities in the 
Weyl tensor. 

For homogeneous FLRW, C = (which keeps the shear term ^ = 0),TZ = — 3/2fim(l + z)^ 
and the affine parameter is related to the redshift by 



For Swiss cheese with Schwarzschild inhomogeneities, (Q) reduces to a single linear and 
homogeneous 3'''^-order differential equation for A, provided the light beam interacts with 
enough randomly distributed inhomogeneities as it travels from source to observer, jsij . This 
3'' -order equation was an analytic version of 3 and which accounted for the decreased 
mass density in the observing beam and for the lensing of galaxies, except it was derived 
for an exact GR cosmology. In |12] the entire calculation was extended to include A 7^ by 
simply including the Qa term in Q. 



(1 + ^)'/rT^u^T^M(rT^F^^ 



dz 



(2) 
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A more universal application of (0), as recognized by Dyer and Roeder occurs when 
the beam interacts infrequently with inhomogeneities or when the interactions introduce 
negligible shear. In this case ^ ~ and ^J~A. satisfies a linear homogeneous second order 
differential equation: 

n. (3) 



VC4 

Values for IZ can be computed for exact gravity fields such as for FLRW as given above and 
for Swiss cheese where 7^ fluctuates, dropping from the FLRW value to zero in voids. For 



Swiss cheese, 



2l| and 



12j |. where the driving term IZ can be replaced by its average along 



the light path the area equation becomes, 

^ = + (4) 

\/ A ^ Pm 

where pH/Pm is the filling fraction (often denoted by a) and is the transparent fraction 
of the total matter density that appears in the beam. Many refer to this equation as the 
Dyer-Roeder equation, however, to do so does an injustice to Zel'dovich who first created 
and solved the equation, to Sachs whose optical scalar equations it is rigorously derived 
from, and to others who developed and solved the equation prior to Dyer and Roeder's first 
pubhcation Q|. 

Modification of (0) and the value of the Ricci term in Q can also be made to account 
br additional matter contents of the FLRW model such as quintessence, etc., see j29|. jl6|. 



41|, Q and Q- 

It should be understood that there are no exact Swiss cheese type 



models containing pressure in which to verify perturbative assumptions such as the affine 
parameter and the radius-redshift remaining as they are in the corresponding homogeneous 
FLRW model. 

In the next section is converted into its Lame' form and in the following section it is 
solved to give D{z) for some special cases. One new result (jHJ, seems particularly useful, 
that of fiat filled-beam FLRW. 



II. THE PARTIALLY FILLED BEAM EQUATION 

In this section the partially filled beam area-redshift equation for pressure-free FLRW is 
revealed to be a standard Lame equation when the dependent and independent variables 
are appropriately redefined. 
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We first rewrite equation (jH) using © as 



:i + Z)V1 + ^mZ + ^k[{1 + - 1] X 

- VtraZ + f^A 

(3 + z/)(2-z/) 



^(1 + ^) Vl + ^mZ + 1^a[(1 + Z)-^ - 1] 



n^{l + zfJA{z) = Q, (5) 
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where the parameter v (the filling parameter) has been introduced to accommodate standard 
mathematical notation for special functions 

ot viv + 1) 

— = ^ \ < z/ < 2 

Pm 6 

— = a = 1 = . (6j 

Pm Pm 6 

Luminosity distance is related to a particular solution A{z) of this equation by 



+ ..)/||. (7) 

A|o is the observed area of a beam of light that started expanding with inital solid angle 
5VLs from a point source at redshift Zs- It is evaluated by integrating equation © from the 
source 2; = 2;^ to the observer z = with initial data which makes the wave front satisfy 
Euclidean geometry when leaving the source (area=radius^ x solid angle): 



where in FLRW the value of the Hubble parameter at Zs is related to the current value Hq 

at z = by 



H, = Ho{l + z,) ^1 + n^z, + fiA[(l + z,)-2 - 1]. (9) 
In j23| equation (jSJ was transformed into the form of a standard Heun equation by changing 



the independent variable from z to y and the dependent variable from y A{z) to h. 

A. Qo = nrn + ^A7^ 1, The Non-Flat FLRW Models 

For the non-flat cases the required change of variables is 
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The resulting equation is 



£h [l + ly)y dh \u{u + l)y+l 



h = 0, (11) 



dy^ + y'^ - ba dy y^ + y"^ - 
where the constant 6^ depends on the two FLRW density parameters Qm and Q\, 

bn = -nl^A/il - ^of. (12) 

Equation (fTT| was recognized as a Heun [3] equation, (see j3| for details) and its solution 
was given in terms of linear Heun functions, HI by Kantowski 23|. Sereno et al. 41 1 
found that it remains Heun if the cosmological constant is replaced by domain wall sources; 
however, when replaced by string networks it reduces to hypergeometric. Equations for 
other types of quintessence sources are yet to be classified. 

Even though the Heun equation is only slightly more complicated than the hypergeometric 
equation (e.g., it has 4 regular singular points rather than 3) Heun functions are not yet 
available in standard computer libraries. Consequently, such expressions are not yet useful 
for comparison with data. Because the exponents of three of the singular points of the area 
equation (in standard Heun form) are and 1/2 [see Eqn. (13) in 2^], equation ()11|) is in fact 
equivalent to the Lame equation, as was pointed out at that time. Details of transforming 
(jllll . i.e., © into Lame form will now be given. The required change of independent variable 
is: 

uiz) = / , < 0, 13 

/4t3 _ t/12 - g. 



where integration is along the positive real axis and 



piz) = , ^"-^ , (1 + ^) - - >--. (14) 
' 4 1 - fin ^ 12 12 ^ ^ 



The constant is given by 



^3 = ^(1-^), (15) 
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where the constant b is related to bn of ()12|) as 



in 0, 



97 02 O 



and the parameter n is determined by the sign of the 3-curvature 



The three cubics appearing in (0), (fTT|) . and (fT^ are related by: 



1 



1 + zfil + n„,z) + nA[i - (1 + z)']}, (18) 



4 / |1 -^^oP 

and are all positive for Big Bang cosmologies. The resulting Lame equation^ in Weierstrass 



form is: 



K i - 



= 0, (19) 



where V {u; g2, gs) is the Weierstrass elliptic function (see Figs. 1 and 2) defined by 

« = / , . 20 

Joo V4t^ - g2t - ^3 

The discriminant in the general case is 

^^gl-TJgl (21) 

which reduces to 

A - (22) 
^" (12)3 ' \^^) 

for and is positive only for < 6 < 2. For these special models the cubic in (pn|) has 

three real roots and for Big Bang cosmologies these are the recollapsing models. For other 

values of b the cubic has only one real root and a pair of complex conjugates ones. 

The Riemann surface for the integrand of ()20|) is double layered because of the square 

root and has three branch points because of the cubic. The two layers of the square root 

make V{u) = V{—u) and the branch points make V{u) doubly periodic. The integral is path 



^ See section 15.2 of |3- Equation p9|) above is not the most general Lame equation. In general it contains 
two additional arbitrary parameters: H rather than k(1 — v{v + 1)/12), and 52 as in H20|l rather than 
1/12. 
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dependent and by circumscribing two branch points the value of the integral is increased by 
a fixed period. By circumscribing another pair of branch points a second period is obtained. 
A third choice of branch points gives a combination of the first two periods. The domain 
for V{u) is usually given as a fundamental parallogram (cell) centered at the origin of the 
complex u plane and generated by two half periods u and u', see Fig. 1. What is needed 
to determine Luminosity-distance from ()19p are values for u in this cell on the negative real 
axis, as determined by (fT^ . In general, V has only one singularity in a fundamental cell, 
a second order pole at m = with vanishing residue, for which V{u) — l/u"^ is not only 
analytic but vanishes at m = 0. To use (fTT?|) to determine luminosity distance the boundary 
conditions (jH} on the area A must be translate to boundary conditions on h: 



h\s = 0, 
dh 
du 



(23) 



i^oy |1 — ^o| 

The resulting luminosity distance is then given by 

D,{zs) = {l + z,)h{um- (24) 
In Section UTTl equation (fTI?|) is integrated for some special cases. 



B. nQ = VLm + ^A = 1, Flat FLRW Models 

For the fiat case the required change of variables for (0) is 

= (1 + z). 



A = + (25) 

and the resulting equation in Heun form is 

dy^ ^ 1/3 + ^2^^^ dy y^ + n^n^ ' ^ ^ 

This equation was described in |23i] but was never actually written down. However, solutions 
to ()26|). linear Heun functions, were used to give an expression for luminosity distance for 



fiat FLRW. When fin = 1, Kantowski and Thomas 



25[ showed that (jSj) was equivalent 
to the Legendre equation and De{z) was given in terms of associated Legendre functions. 
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as well as equivalent hypergeometric functions. An equivalent 



h ype rgeometric result was 



42 1 have shown that the 



independently given by Demianski et al. Sereno et al. 
equation remains hypergeometric even when generic quintessence is introduced. 

To transform (j26j) into Lame' form the change of dependent variable is the same as in 
(I25|) : however, the change of independent variable from z to u simplifies to 

/^(^^ dt 1 ^ /I 1 7 1 \ 

' VWT4: Jp(z) \6 2 6 p{z^^ ' 



oo 



where 

piz) ^ + z) >0. (28) 

The area equation for fiat FLRW in Lame' form becomes: 

"^'^ u{u + l)V{u;0,-A)h = 0, (29) 



where 



du 

Boundary conditions on h change to: 



du v/TO^-^Vl + ^U(l + ^)^-l] (30) 



h\s = 0, 
dh 



du 



2c 



(31) 



but the connection of h to Dii{z) is unchanged from 



III. THREE SOLVABLE CASES 

To give useful and/or interesting expressions for luminosity distance as a function of 
redshift requires a knowledge of the special functions associated with the defining differential 
equations (fTTj) or (fT^. In Kantowski [2^ special functions associated with (fTT|) called Heun 
functions were used to give Df {z); however, Lame' functions defined by the Lame' equation 
(|19|) are not defined for arbitrary values of the constant H (see footnote 1). Most of the 
special functions defined by (jl9p relate to boundary value problems where H is restricted to 
discrete characteristic values, see Q, Q, and Q The characteristic values do not coincide 
with the values needed in p9|) . However, there are three values {0, 1,2} for u for which a 
complete solution to fT^ can be given. InQ it was recognized that the solution of (0) could 
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be given as elliptic integrals in Legendre form for these three cases. It is for these same three 
cases that the complete solution to (fT^ can be given and hence for these three cases the 
luminosity distance can be given by additionally applying boundary conditions, or (jSH), 
and using (j24j) . Because the Weierstrass form of Lame equation contains the Weierstrass 
elliptic function its solution will be given in terms of Weierstrass elliptic integrals. 



A. 1/ = 0, The Filled Beam Case, Standard FLRW 

The one case for which (jl9j) is easily solvable is z/ = 0, the homogeneous filled beam 
FLRW case. 



1. = 0, ilo = + 7^ 1) The Non-Flat Case 



The general solution is: 

h{u) = ci S^,[u - Co], (32) 

where Cq and Ci are integration constants and S^[] is one of two functions, 

sinh[ ] : K, = —1, 
sin[ ] : K = +1. 
Eqn. (|32|l results in an expression for luminosity distance 

D^{Q^, fiA, u = 0;z)= [u{z) - n(0)] , (33) 



Hoy\l — 

where the source redshift has been simplified from Zg to z and u{z) is defined by p3|) . If 

V^^ (p ; 92, ds), the inverse of the Weierstrass function defined in 

dt 



v~'{p; 92,93) = r 

J oo 



V4t3 - g2t - gs 

is used to evaluate u{z) the two inverse Weierstrass functions in can be combined using 
the identity: 



V(u + v) 



V'{u) - V'{v) 
V{u) - V{v) 



V{u)-V{v), 



to give 

u{z) -m(0) 
/ 



V-'[v{zY^,93)-V-'[p{Q);^,93 



(35) 



(36) 



-V 



-1 



^(1 + Zfil + VlraZ) - Z{Z + 2)Vlf, + 1 


' - n^{z + 2)z^ 


4|1 - fio 


z^ 



■^^^9z I > 0. 
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2611 and 



27|.2 Care 



The resulting expression for ()33|) was first given by Kaufman Schucking 
has to be exercised when using the identity (jHBj) in ()33|) . The inverse function ()34|) is multi- 
valued (e.g., its value (±) depends on p 's location on the Riemann surface of the square root 
appearing in (j34j) ) . In 2j|, (j33|) was given in terms of incomplete Legendre elliptic integral 
of the first kind, F(0,k), rather than inverse Weierstrass elliptic integral V^^{p; g2, gs)- 



2. = 0, flo = + = 1) Standard FLRW, The Filled-Beam Flat Case 



For the fiat case the general solution to ()29|) is 



h{u) =ci{u- Co), 



which results in a luminosity-distance 



2c(l + z) 



(37) 



(38) 



The z dependence of u{z) — u{0) can be evaluated using which is equivalent to a 
hypergeometric function by ^T7\. 

u{z) = v-'(p{zy,o,-4 

1 ^ fl 1,7 1 



6 2 6 p{z) 



(39) 



where p{z) is given in ()28p. The third equality is simply a relation that exists between all 
associated Legendre functions and some hypergeometric functions. The addition law (j35| 
results in the simplification 



u{z) - m(0) = {p{z)- 0, -4) - (p(0); 0, -4) 



2Jl + Vt^z{^ -I- 3^ + z2) + 2 + Vt^z(?, + z) 



;0,-4 > 0.(40) 



The expression obtained by using the last equality in pUj) results in an expression obtained by 
when the identity V {u; g2, gz)=T'^V (Tu; T'^g2, T~^g^) is used. When the equivalence 



^ A factor of 1/2 was left out of (4) of compared to (32) of [2^ and compared to (|33|l combined with 
above. 
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of the inverse Weierstrass function V~^{p;0, —4) and the hypergeometric function 2-^1 as 
imphed by is used in (jHHj) a new and simphfied expression for luminosity distance for 
standard (filled beam) flat FLRW results: 

2c(l + z)z 



1 - Vtm, V 



0;^ 



1 17 

6'2' 6' 



2 + fi^z(3 + z) 

3 



1/2 



X 



2^1 + fi„2(3 + 3z + + 2 + fi^2(3 + z) 

c 2^/6 r (I) (1 + 2) 2 
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z2) + 2 + r]^z(3 + ^) 



1/4 



X 



2JI + 1]„2(3 + 32 + + 2 + fi„2(3 + z) 



3\ 



(41) 



25 1, and jol for previously known expressions. The first expression in ()41|) should 
be quite efficient for numerical evaluation because 2-^1 = 1 + 0{z^) and the argument of 
the hypergeometric function remains close to zero for reasonable z. The second expres- 
sion involving the associated Legendre function -PJ^/g seem less useful but is included for 
completeness. 



B. iy=l, The 66% Filled Beam FLRW 



This is a partially filled beam case corresponding to having i/(l + z/)/6 = 1/3 of the mass 
density excluded from the observing beam. For this case (jl9|) can be integrated in a straight 
forward manner giving the general solution in terms of two integration constants Cq and Ci'. 

du 



\ 



486 - h 



3325 



CO 



^3 



(42) 



where 



sinh[ ] : 6 < 
] = < sin[ ] : < 6 < 486, 
sinh[ 1 : 486 < h. 

n 

Special cases like b = 486 were dealt with in [2^ but are simply neglected here. The above 
result can be checked using the identities 



V" = 



:92, 
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{V'f = _ g^v - 93 



(43) 



satisfied by all by all Weierstrass V functions. The denominator satisfies: 

+ V{u-—,g3) > — (44) 



6 V 12 

and can cause complications for negative curvature models if it vanishes along the path of 
integration. When boundary conditions are applied, the luminosity distance becomes: 



Ho 
Sign 



l + z] 



\ 



[3 - fi„(l + z)/{l - Qo)] [3 - nj{l - Qo)] 



f+piO) 



:i-r]o)(486 -6)/(332) 



X 



(Slm,fiA>z) 



486 -6 



3325 



du 



n(0) 



where u{z) is given by (fT^ and 



5, 



cosh[P ] : 6 < & [3 - Qm{l + z)/{l - Qo)] [3 - a^/(l - Qo)] < 0, 

sinh[ ] : 6 < & [3 - ^]™(l + z)/{l - l^o)] [3 - ^^m/(l - ^^o)] > 0, 

sin[ ] : 0<b< 486, 

sinhf 1 : 486 < b. 



For the first case, 6 < 0, the principal value of the integral, P/, in ()45|1 is needed when 
the denominator of the integrand vanishes between the source and observer, see (0H). In 
this case the argument of hyperbolic trig function sinh[] develops an imaginary part, zvr/2, 
which turns it into a cosh[] of the principal part. 

The integral in is an elliptic integral written in terms of the Weierstrass V function 
and as such can always be evaluated in terms of the three Weierstrass functions V~^, (, and 
a (see 13.14 of ^^). For ()45j) the integrand is an elliptic function with the same fundamental 
cell as V{u; 1/12, (73) and has two first order poles at u = j3±, {(3- = — /3+) within the cell 
shown in Figs. 1 and 2. The points P± are those two lattice points defined by 



P = V-\-K5/6;l/12,g3) 



(46) 



and which lie within the fundamental cell. An elliptic function is uniquely given (up to a 
constant that can easily be evaluated) by the principal parts of its Laurent expansions about 
its poles within a fundamental cell. Because the principal parts can be expressed in terms 
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of the Weierstrass ^-function and its derivatives, a unique expansion of the elhptic function 
exists in terms of these functions. The expansion for the integrand in (pSjl is: 
1 1 



V 



12 



2 C gsj +c(u- p+; ^, - C + P+; ^3) 

(47) 

In fl47|) the first two terms are constants and the second two are the principal parts at the 
two first order poles in the fundamental cell. The value of 

1 



V f3. 



12 



,93 



-k(486 - b) 
63 



(48) 



is real except when < 6 < 486. When integrated (jlTjl gives 



u{z) 
u{0) 



du 



+ log 



a {u{z) 






a (m(0) 






a (m(0) 
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a {u{z) 







(49) 

The simplification of ()49j) similar to the simplification of ()56|1 must exist but escapes the 
author. To evaluate and hence (jlH]) the three quantities 'u(z),-u(0), and u{z) — 'u(O) can 
be evaluated as in ()36|) combined with (fT^ . Fortran routines for C, and a are not available 



at this time, however, they do exist in Mathematica. In |2J| the area equation (jSJ was 
integrated directly giving luminosity-distance in terms of Legendre elliptic integrals -F(0, k) 
and n(0, a^, k). Fortran routines for these are available (see the Numerical Recipes website). 



C. V = 2, The Empty Beam Case 



This is the case where all material is exterior to the observation beam and its effect 
through the Weyl curvature (lensing) is neglected. For this case (fT^ can be integrated in a 
straight forward manner giving the general solution in terms of two integration constants Cq 
and Ci, 



h 



Cl 



1^ 1 
12 V 12'^^ 



du 



CO 



2 • 



(50) 



'f^+v{u-l^,g,] 

As with the u = 1 case, this result can be checked by using ()4H|1 . When boundary conditions 
(I2H|1 are applied the luminosity distance becomes: 



De{Qrn,^A,i' = 2; 2;) 



c 



411 -^]n 



+ 



1-fin 



u{z) 
«{0) 



du 



f^+Viu;^,g, 



12' 



2 ' 



(51) 
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where u{z) is given by ()13p. For Big Bang models the denominator in the integral doesn't 
vanish for any z > and hence problems that occurred for the u = 1 case do not appear for 
u = 2. Just as with (|45p the integral in (|5Up is an elliptic integral written in terms of the 
Weierstrass V functions and as such can also be evaluated in terms of the three Weierstrass 
function V~^,(j ^ind a. For (|51|) the integrand is again an elliptic function with the same 
fundamental cell as V{u; 1/12, g^) but has two second order poles at m = P±, (/?_ = — /?+) 
within this cell. The two poles /3± now belong to the lattice defined by 



(3 = V-\-K/12;l/12,g3). 



(52) 



The unique expansion for the integrand in (j^Tj) is: 
1 1 



[C'(-/3+; 1/12, ^73) + C'(/3+; 1/12,^73)] 



- [C' {u - I3+- 1/12, g^) + C'{u + (3+- 1/12, (73)] }, 



(53) 



where C'l"^) = is the first derivative of the Weierstrass (^-function. In (j53p the first 

two terms are constants and the second two are the principal parts at the two second order 
poles. 

This expression can be easily integrated to give a value for (jS^, 



c \/\l —Qo\ ,^ r 

C (u{z) - ^,g,^ - C (ix(O) - g, 
C (u{z) + ^, (73) - C (w(0) + P+; ^, gs 



(54) 



The ( terms can be combined using the ^-identity: 



1 



({u + v) = C{u) + C{v) + - 



V'iu) - V'iy) 
V{u) - V{y) 



(55) 



along with Using P(/?+) = -k/12 and P'(/?+) = ±^Khl{?,n^) = ±fi„yni/(4|l 

f2o|^''^) the resulting luminosity-distance becomes 



Di{Qrn,^A,T^ = '2] z) 



i-no\ 



K 

+ X + 



1 



2 4|l-fio| 



■ / 1 JIh 
2( — + ^ 



^z^ {1 + zy 



[l + zY ^ [u{z)-um-2c(u{z)-u{0y,^,g3 

^{i + zyii + nmz)-z{z + 2)n^ 
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+ 3 



+ < 



K 



z Z^ 

1 



z (1 + z) 



1 



1/2 



411 - U 



0| 



2;^ 



:i + ^) 



1 + 2)2(1 + 

2 2 \ 



2(2 + 2)11a 

1/2- 



(56) 



The n(z) — m(0) arguments are given in (jHBj) resulting in a distance-redshift expressed in 
terms of the two independent Weierstrass elhptic integrals and 

If the integration in (j5T|l is changed from m to 2; by using (|T3|l and (|THl . the luminosity 
distance becomes the familiar elliptic integral, 

dz 



1 + z) 



1 + Z) V(l + + ^mZ) - Z{Z + 2)fi^ 



(57) 



obtainable directly from (jSj). In 2J] the value of this integral was given in terms of Legendre 
elliptic integrals F(0,k) and E{(f),k) in a somewhat more complicated expression. 



IV. CONCLUSIONS 

In this paper the partially filled beam Hubble relation for pressure free FLRW cosmologies 
has been shown to satisfy the Lame' equation ()19|) or ()29p when appropriate dependent 
and independent variables are chosen. For three values of the beam's filling parameter 
{u = 0, 1,2), general solutions to the Lame' equation were given in terms of Weierstrass 
elliptic integrals. These integrals were used to give the Hubble relation for these three 
cases. Practical use of these expressions awaits implementation of fast fortran routines for 
the Weierstrass integrals. The results given here were checked using the slow routines of 
Mathematica. The Lame' equation has not found its way into the physics or astron om^y 
literature much in the past, however, recently some use in string theory is showing up. 



Lomy 



The partially filled 
Legendre equation 



J3eam Hubble relation for flat FLRW is known to satisfy the associated 
and hence can be written as a linear combination of two associated 
Legendre functions. In general the two coefficients of these functions both depend on asso- 
ciated Legendre functions bringing the total to four. For the filled beam case, two of the 
functions reduce to ordinary functions leaving two to be evaluated as associated Legendre or 
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as their equivalent hypergeometric functions. For this special case an addition formula for 
the Weierstrass V functions allowed the discovery of a Hubble relation (PT|) containing only 
one associated Legendre or equivalent hypergeometric function. A fortran implementation 



of this expression should rapidly converge (see [http:// www. nhn . ou . edu / ~thomas /z2dl . html 
for earlier code). 

The Hubble relations given here are good for observations made on distant objects when 
lensing by inhomogeneous mass external to the observing beam can be neglected and when 
transparent mass interior to the beam is scattered along the beam and not concentrated at 
isolated regions. For cases where mass concentrations external and/or internal are impor- 
tant, expressions given here are necessary to compute source-deflector, source-observer, etc., 
distances for gravitational lensing, see jj] and . In 2^ it was argued that even though the 
filling parameter v could modify observed magnitudes by small amounts, ~ 0.05, a tenfold 
increase the number of SNela m{z) data (~ 600) would statistically necessitate limiting its 
value. This amount of data is to become available from the SuperNova Acceleration Probe 
(SNAP - p^p://snap.lbl.gov |. 



The miz) relations given here can also be used to help understand ray tracing simulations 

n n n 

;hat predict weak lensing amplifications distributions for SNe, e.g., j46||, j36|], [18||, and 
^. The minimum amplification for weak lensing occurs for beams that have a minimum 
amount of transparent and homogeneous mass within them and for which inhomogeneities 
are far enough away to have negligible lensing effect (see Fig. 3 for a generic histogram of 
lensing amplifications, f47'l). Some work has been done to find an effective filling fraction 

n n 

a = Ph/ Pm = (3 + z^)(2 — z/)/6, |4^- Such a filling parameter Veff is not necessarily 
available a priori for the simulation; however, it can be obtained by comparison of the 
simulation's histogram with the partially filled beam de-amplifications given here, 

Amp^i^ = Dj^rim, ^A, iy = 0; z)/Dj{Vt^, ^a, v^ff] z). (58) 

The ratio of the luminosity with an effective filling parameter, to the average luminosity can 
be taken to be the minimum amplification of the simulation. 
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FIG. 1: A typical fundamental cell (domain) for the doubly periodic Weierstrass elliptic function 
7^ (-u; 1/12,53 = -1.75). The cell and half periods, u = -1.39i and u' = 1.21 + 0.69z, are for a 
universe with = 0.2 and CIa = 0.7. The periodic locations of the 2"'^-order poles are indicated 
by the X's. 
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FIG. 2: Plots of V {u; 1/12, = —1.75) along the real and imaginary u-axies for a universe with 
Um = 0.2 and Q,a = 0.7. The function is real on these lines. 
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FIG. 3: A schematic distribution of weak lensing amplifications (relative to a homogeneous uni- 
verse) expected for observing a set of SNe at a fixed redshift. The indicated minimum can be used 
to determine an effective value of the filling parameter v by using H58() . 
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